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ABSTRACT 

Rotation curves of dark-mattcr-dominated galaxies measure the mass profiles 
of galactic halos and thereby test theories of their cosmological origin. While 
attention has focused lately on the possible discrepancy at small galactocentric 
radii between observed rotation curves and the singular density profiles predicted 
by N-body simulations of the Cold Dark Matter (CDM) model, the observed ro- 
tation curves nevertheless contain valuable additional information with which to 
test the theory and constrain the fundamental cosmological parameters, despite 
this uncertainty at small radii. An analytical model we derived elsewhere for the 
postcoUapse equilibrium of cosmological halos as truncated, nonsingular, isother- 
mal spheres (TIS) reproduces many of the average properties of halos in CDM 
simulations to good accuracy, including the density profiles outside the central 
region. The circular velocity profile of this TIS model is, moreover, in excellent 
agreement with the observed ones and yields the mass and formation epoch of an 
observed halo from the parameters of its rotation curve. This allows us to predict 
correlations amongst rotation curve parameters, such as the maximum velocity 
and the radius at which it occurs, for different mass halos forming at different 
epochs in the CDM model. As an example, we derive the observed I'max-'^max 
relation analytically, with preference for the flat ACDM model. 

Subject headings: cosmology: theory - dark matter - galaxies: general - 
galaxies: formation - galaxies: halos - galaxies: kinematics and dynamics 
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1. Introduction 

The rotation curves of dark-matter-dominated galaxies probe the mass profile of their 
galactic halos. The observed rotation curves of dwarf and low-surface-brightness (LSB) disk 
galaxies, therefore, offer a relatively direct test of theories of cosmological halo formation, 
free of the complicating effects of the dynamical coupling of dissipationless dark matter and 
dissipative baryonic matter which affect the mass profiles of baryon-dominated galaxies. 
Recently, attention has focused on the apparent conflict between the singular halo density 
profiles predicted by CDM N-body simulations and the observed rotation curves, which 
favor a flat-density core (cf. Moore et al. 1999). This has led to intense scrutiny of both 
the observations and the CDM model, primarily in one of three directions: improvement 
in the numerical resolving power of the CDM N-body simulations to determine better the 
logarithmic slope of the predicted density profiles at small radii (e.g. Moore et al. 1999), 
ideas for modifying the microscopic properties of the dark matter so as to retain the more 
successful aspects of the CDM model while flattening the halo density profiles at small radii 
(e.g. Dave et al. 2000, and references therein), and suggestions that the rotation curve data 
lack sufficient spatial resolution near the center to distinguish unambiguously between a 
density profile with a fiat-density core and the singular profiles predicted by CDM N-body 
simulations (e.g. van den Bosch & Swaters 2000). Despite these uncertainties, a meaningful 
comparison between the data and predictions for the global properties of the rotation curves 
of dark-matter-dominated galaxies is possible which serves to test the CDM model and 
discriminate amongst different background cosmologies. To demonstrate this, wc show how 
the mass and formation epoch of a halo can be extracted from its observed rotation curve, 
by application of an analytical model for halo formation which reproduces the empirical 
description of these rotation curves extremely well. 
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Burkert (1995) showed that the observed rotation curves of several dark-matter- 
dominated dwarf galaxies are consistent with a common density profile with a flat-density 
core, according to 

(r/ro,B + l)(rV4B + l)' ^ ' 

Kravtsov et al. (1998) found that the rotation curves for a larger sample which included 

both dwarf and LSB galaxies are well fit by a similar universal halo profile given by 

P^""^ = (r/r5)^(l + (r/r5)")(/3-^)/°' ^'^^ 
with (a,/9, 7) = (2,3,0.2), which is only slightly steeper than the Burkert profile in the 
core but shares the slope at large radii, where p oc . The circular velocity profiles [i.e. 
v{r) = {GM[< r]/r)^/^] for equations (1) and (2) are virtually indistinguishable at all radii, 
even near the center, where v oc r^"'"'/^, since 1 — 7/2 = 1 or 0.9, respectively. 

This contrasts with the halo profiles found by N-body simulations of the standard 
CDM model. The universal fitting formula for simulated halos reported by Navarro, Prenk, 
& White (1997; NFW) is equation (2) with (a,/3,7) = (1,3, 1), while Moore et al. (1999) 
find (a,/?, 7) = (1.5,3,1.5). However, the circular velocity profiles for these halos differ 
significantly from those implied by equations (1) and (2) only at very small radii where the 
uncertainties in the shape of the observed rotation curves currently makes discrimination 
difficult. In the meantime, equation (1) generally fits the data better than does the NFW 
profile even when the latter is also an acceptable fit. 

The Burkert profile, then, continues to serve as a useful empirical description of the 
universal mass profiles of dark-matter-dominated galactic halos. We show here that the 
truncated, nonsingular isothermal sphere model we derived elsewhere (Shapiro, Iliev & 
Raga 1998, Paper 1; Iliev & Shapiro 2000, Paper 11) has a density profile for which the 
circular velocity profile is essentially indistinguishable from that of the Burkert profile and, 
as such, provides a theoretical motivation for the latter. The TIS model goes well beyond 
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the prediction of density profile and rotation curve, however, to provide the size, mass, 
velocity dispersion and collapse epoch of the halo, as well. This makes possible the further 
interpretation of observations of dark-matter-dominated galaxies for comparison with the 
predictions of various cosmological models. At the same time, the TIS model can be 
shown to reproduce many of the average properties of the halos found in simulations of the 
standard CDM model, outside of the innermost region where the TIS halo has a fiat-density 
core, unlike the CDM halos. As a result, a comparison of the analytical TIS predictions with 
observed properties of galaxies also offers insight into the standard CDM model. Finally, 
if suggestions hke the self-interacting dark matter proposal of Spergel & Steinhardt (2000) 
are correct, that CDM might be more "coUisional" as a way to eliminate the central cusp of 
standard CDM halos, then our TIS solution will also apply to these models, to the extent 
that the halo relaxation process makes the final equilibrium approximately isothermal. 

In § 2, we briefly summarize the relevant properties of the TIS model and present 
a simple analytical formula for it with which to fit an observed rotation curve. In § 3 
we demonstrate the excellent agreement between the TIS model rotation curve and that 
implied by equation (1) and show how this allows us to deduce the total mass and collapse 
epoch of a given halo directly from the parameters of its observed rotation curve. In § 4, 
we describe how the dependence of the average formation epoch of a halo on its mass in 
the CDM model results in statistical correlations amongst the parameters of the observed 
rotation curves. As an example, we use this approach to derive analytically the known 
correlation between the maximum velocity of each rotation curve and the radius at which 
it occurs. 
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2. The Truncated Isothermal Sphere Model 

The TIS model is a particular solution of the Lane-Emden equation (suitably modified 
when A 7^ 0) which results from the collapse and virialization of a top-hat density 
perturbation (c.f. Paper I for Einstein-de Sitter [EdS] universe and Paper II for Qq < 1 
and Ao = or 1 — Qo)- The size rt and velocity dispersion ay are unique functions of the 
mass and redshift of formation of the object for a given background universe. While the 
Lane-Emden equation requires a straightforward numerical solution, Paper I provides a 
convenient analytical fitting formula, 

A B 



P{r) = Po 



(3) 



where C = r/ro and {A, a^ B, 6^) = (21.38, 9.08, 19.81, 14.62), accurate to within 3% over 
the full range < C < Ct ~ 30 for both the EdS case and most low-density models of 
interest (i.e. flo > 0.3). [Note: Our definition of the core radius is ro,Tis = ?^King/3, where 
TKing is the "King radius" defined in Binney & Tremaine (1987), p. 228.] Equation (3) can 
be integrated to yield an analytical fitting formula for the TIS rotation curve, as well, given 
by 



^ ^ ' A-B + 



ftStan"^ ( ^1 -aAtan"^ ' ^ 



1/2 



(4) 



(^v [ C 

where ay — (47rG'po^o)^^^- This fit has a fractional error less than 1% over the full range of 
radii < r < for all matter-dominated background models. With a nonzero cosmological 
constant, the circular velocity becomes v{r) — {GM[< r]/r)^/^[l — 2px/p{r)Y^'^, where px is 
the constant vacuum energy density associated with the cosmological constant. For cases 
of current interest for a fiat universe with A 7^ (i.e. Qq > 0.3), equation (4) can still be 
used, however, since it departs from the exact solution only slightly in the outer halo (i.e. 
at r ~ rt) for halos collapsing even as late as 2; = 0; the fractional error is less than 6% at 
all radn and less than 1% for r < (2/3 )rt. 

The TIS model quantitatively reproduces the average structural properties of halos 
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found in CDM simulations to good accuracy, suggesting that it is a useful analytical 
approximation for halos which form from more realistic initial conditions. An exception to 
this agreement is the very inner profile, where the TIS has a uniform- density core instead 
of a central cusp. Our TIS predictions agree to astonishingly high accuracy (i.e. to of 
order 1%) with the cluster mass-radius and radius-temperature relationships and integrated 
mass profiles derived from detailed CDM simulations of X-ray cluster formation by Evrard, 
Metzler, and Navarro (1996). Apparently, these simulation results are not sensitive to our 
disagreement in the core. A direct comparison of the TIS and NFW mass profiles reveals 
a very close agreement (fractional deviation of less than ~ 10%) at all radii outside of 
a few TIS core radii (i.e. about one King radius), for NFW concentration parameters 
4 < cnfw ^ 7. 

3. Application to Galaxy Rotation Curves 

The rotation curve for equation (1) is given by 

VB{r) _ f ln[(CB + l)^(C| + l)]-2tan-^(CB) 

'V*,B I Cb 

where (b = r/ro^B and v^^b = (ttGpo.b^^o.b)^^^- To compare our TIS directly with 
equation (1), we solve for the ratios Po,tis/Po,b and ro,Tis/'"o,B which minimize the oi 
the fit of equation (5) to equation (4) over radii from to r^. The result in Figure 1, 
with po,Tis/po,B = 0.7790 and ro,TisAo,B — 0.3276, shows extremely good agreement, 
with a fractional deviation below 10% from O.Olr^ to rt and below 4% over the range 
r > 0.03rt ~ ro,Tis- For the TIS, the maximum circular velocity and its location are 
'i'max,Tis = 1.5867(Tv,Tis ^t Tmax.Tis = 8.99ro,Tis for all matter-dominated cosmologies. For a 
flat universe with A 7^ 0, {Qq > 0.3), these numbers depend weakly upon collapse redshift, 
but are reduced by no more than 0.2% and 1.7%, respectively, for Zcou = 0, and by even 
less for earher collapse. For the Burkert profile, t'max.B = 1.2143i'*,b at Tmax.B = 3.2446 ro,B- 
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Hence, our best fit finds that Vmax and rmax for the TIS and Burkert profiles are extremely 

close, with ^^max.TIs/^'max.B = 1-02 and rmax,TIs/^^max,B = 0.91. 

In short, our TIS model provides a solid, theoretical underpinning for the empirical 
fitting formula of Burkert, and, by extension, a self-consistent theoretical explanation of the 
observed galaxy rotation curves it was invented to fit. In addition, the TIS model allows us 
to calculate the total mass Mq, collapse epoch ^cou, and other parameters of each observed 
dark-matter-dominated halo from its rotation curve, as follows. Assuming Vmax and Tmax 
are provided by observation^, the mass of the galaxy can be shown to be 



for any matter-dominated universe. When A 7^ 0, Mq oc rmax^max; but the coefficient 
depends on ^icou and background cosmology (Paper II). For any flat universe (A 7^ 0) of 
current interest, equation (6) is still a very good approximation. For D,o — 1 — Xq > 0.3, 
equation (6) underestimates the mass by less than 6.3%, 2.8%, and 1.2% for Zcow — 0,0.5, 
and 1, respectively. 

In terms of Tmax and I'max, -^cou is given to good accuracy by the implicit equation 



where F = {[Qo/Q(^coii)][Ac,sus/187r2]}V3(i + ^^^^)^ q(^) ^ + ^)3]/[(i _ _ x^){i + 



zY + Qo(l + + Ao], and Ac,sus is the density [in units of Pcrit(-2coii)] after top-hat collapse 
and virialization for the standard uniform sphere approximation (SUS)^. For the EdS case, 

^In practice, observations are generally restricted to the inner parts of galaxy rotation 
curves, often not extending to radii as large as rmax- In that case, I'max and rmax are inferred 
by fitting the observed rotation curve to the TIS rotation curve at other radii. 

^This Ac,sus is well-approximated by Ac,sus = IStt^ -|- cix — C2x'^, where x = ^(zcoii) — 1, 
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F — {1 + 2;coii), while for open, matter-dominated and flat cases, F — > Qq ' (1 + ^^cou) at 
early times [i.e. x ^ 0]. The coefficient on the r.h.s. of equation (7) is correct for any 
matter- dominated cosmology and a very good approximation for flat (A ^ 0) cases of 
current interest. For Qq > 0.3, it underestimates (1 -|- Zcow) by less than 2.5% for 2;coii > 0. 

4. Statistical Correlations Amongst the Properties of Dwarf and LSB 

Galaixies: The fmax-^'max Relation 

The halos in our TIS model are fully described for a given background universe by 
their total mass Mq and collapse epoch ^coii- In hierarchical models of structure formation 
like CDM, these are not statistically independent parameters, however. Smaller mass 
halos on average collapse earlier and are denser than larger mass halos. In terms of 
galaxy rotation-curve parameters, this dependence should be observable, for example, as 
a correlation between Vj^ax and Tmax- Mori & Burkert (2000) used the Burkert (1995) flts 
to rotation curves of dwarf galaxies to report such an observed correlation, expressed as 
follows: 

^^max.B = 9.81 (rmax,B/lkpc)^^^ kms"^ (8) 

Kravtsov et al. (1998) also found a correlation using flts to observed rotation curves. Their 
results for a sample of dwarf and LSB galaxies are shown in Figure 2, along with that in 
equation (8). Kravtsov et al. (1998) further showed that the results of their CDM N-body 
simulations agreed with these data points.^ This suggests that if we can relate the mass and 

and ci = 82 (60) and C2 = 39 (32) for the flat (open) cases, Qq + — ^ (^o < 1, Aq = 0), 
respectively (Bryan & Norman 1998). 

^This result by Kravtsov et al. (1998) is not sensitive to the question of whether their 
simulations adequately resolved the halo density proflles at very small radii. 
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coUapse epoch of our TIS model halos in a statistical way within the context of the CDM 
model, a comparison of our predicted Vmax-^max relation with this observed one will further 
check the relevance of our TIS model to halo formation from realistic initial conditions 
and, at the same time, give a theoretical explanation for both the data and the simulation 
results. 

To predict the Vmax-rmax correlation for a CDM universe using our TIS model, we apply 
the well-known Press-Schechter (PS) approximation to derive Zcon{Mo) - the typical collapse 
epoch for a halo with a given mass. Halos of mass M which collapse when o"(M) = ^crit/'^ 
are referred to as 'V-cr" fluctuations, where a{M) is the standard deviation of the density 
fluctuations at Zcou according to linear perturbation theory, after the density field is filtered 
on the scale M, and ^crit is the amplitude of a top-hat perturbation according to linear 
theory at the epoch Zcoii at which the exact solution predicts infinite collapse. The typical 
collapse epoch for halos of a given mass is that for which u — I, the 1-a fluctuations. 

Our results for 1-a fluctuations are shown in Figure 2 (upper panel), for different 
background cosmologies. The flat, untilted and the open, slightly-tilted {up — 1.14) models 
are in reasonable agreement with the observed I'max — '"max relation, while the untilted and 
strongly tilted {up = 1.3) open models are not. The empirical Burkert scaling relation is 
closely approximated by the currently-favored ACDM model, less well by other models. 

In practice, observed galaxies should exhibit a statistical spread of halo properties in 
accord with the expectations of the Gaussian statistics of the density fluctuations which 
formed them. Since observed galaxies will not all be "typical" , the scatter of the data points 
in Figure 2 is natural. To probe this in our model, we calculate the masses and collapse 
redshifts for halos formed by v-a fluctuations for different values of v. As shown in Figure 2 
(lower panel) for ACDM, the Burkert scaling relation is closest to the TIS model prediction 
for 1-a fluctuations, while all the observed galaxy data points except one correspond to 
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0.7 < ^ 1.5. Hence, on average, the galaxies which constitute the observed I'max-'^max 
correlation correspond to halos which formed at close to the typical collapse time expected 
theoretically for objects of that mass. This means that our TIS model is a self- consistent 
explanation for the observed I'max-'^max correlation. 

This success of the TIS model in explaining the observed Vmax-'^max relation and, by 
extension, the CDM simulation results of Kravtsov et al. (1998) which follow it can be 
understood by a completely analytical argument, as follows. We approximate the density 
fluctuation power spectrum as a power-law in wavenumber k, P{k) oc k"'. If we define a 
mass M which corresponds to k according to M oc k~^, then a{M) oc M"*^^"*""^/^ if we set 
n — rieff = — 3(2yefr + 1), where yeff = idlna/d\nM)exaxA a-t the relevant mass scale M. 
Our results for the ACDM case indicate that the galaxies which make up the Vmax — ^max 
data points in Figure 2 collapsed at redshifts 1 < 2;coii ^ 6 with masses in the range 
8 X 10^ < Mo/(Mq/i"^) < 3 X 10^^. Hence, the precollapse fluctuation growth rate is 
approximately EdS, and we can let ^(zcoii) = 1 in equations (6) and (7). In that case, 
(1 + ^cou) OC a{M) oc M-(3+")/6, oc M^^+^'^/^n'^^^ and v^eoc oc M^^-'^^/^^^J/^ which 
combine to yield 

V -V ir Ir ^(i-")/P(5+n)l f9^ 

I'max — '^max,* V' max/ ' max,*^ i V / 

where Vmax,* and rmax,* are for a l-cr fluctuation of flducial mass M*, as given by equations (6) 
and (7) with (1 ^^cou) = (1 + ^rec)cr(-^*, ^rec)/(^crit, whcrc 5crit = 1.6865 and a{M^,Zy.^) 
is the value of (t{M^) evaluated at the epoch of recombination Zrec early enough 
that (1 -|- z)a is independent of z\. Over the relevant mass range M = 10^^^^h~^MQ, 
rieff ^ -2.4 ± 0.1 for our ACDM case. For = lO^o/i-^M©, our COBE-normahzed, flat 
ACDM case (Qo = 0.3, h = 0.65) yields (1 + Zrec)cr{M^, -Zrec) = 5.563, so (1 + ^cou) = 3.30, 
Vmax,* = 53.2 km s~^ and rmax,* = 5.59/i~^kpc. With these values and Ues = —2.4, 
equation (9) yields the TIS model analytical prediction v^ax = (13.0kms~^)(rinax/l kpc)°'^^, 
remarkably close to the empirical relation in equation (8). 
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Fig. 1. — Circular velocity profiles v{r) (in units of ctv^tis) versus radius r (in units of Tq^tis) 
for the TIS model halo (matter-dominated universe) (solid line) and for the best-fit Burkert 
(1995) profile (dotted line), fit over the entire range < r < (i.e. < C < Ct ~ 30). Top 
panel shows the fractional deviation (t'Tis ~ vb)/vtis- 
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Fig. 2. — (a) (upper panel) I'max-'^max correlation. Observed dwarf galaxies (triangles) and 
LSB galaxies (squares) from Kravtsov et al. (1998). Curve labelled "Burkert fit" is equa- 
tion (8) based on a fit by Burkert (1995) to data. TIS+PS predictions are plotted for 
different CDM cosmologies, as labelled: SCDM (Qq = 1, cluster-normalized: (Jgh-i = 0.5), 
OCDM (Qq — 0.3, Ao = 0, COBE-normalized, both untilted [primordial power spectrum 
index Hp = 1] and tilted [Up = 1.14 and 1.3] (subscripts indicate np-value); and ACDM 
(Qo = 1-Ao = 0.3, COBE-normalized). In the untilted ACDM and strongly-tilted {up = 1.3) 
open cases, COBE- and cluster-normahzations are the same. All models assume h — 0.65. 
(b) (lower panel) Same as (a), except TIS-I-PS curves are for the ACDM model only, for fluc- 
tuations of different initial amplitudes for each mass scale, as labelled by the 'V-cr" values: 
l-(7 (i.e. "typical" fluctuations), 0.5-(7,0.7-(7,1.5-(7, and 2-a. 
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